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Abstract 

The reorientation dynamics of local tangent vectors of chains in isotropic 
amorphous melts containing semiflexible model polymers was studied 
by molecular dynamics simulations. The reorientation is strongly in- 
fluenced both by the local chain stiffness and by the overall chain 
length. It takes place by two different subsequent processes: A short- 
time non-exponential decay and a long-time exponential reorientation 
arising from the relaxation of medium-size chain segments. Both pro- 
cesses depend on stiffness and chain length. The strong influence of the 
chain length on the chain dynamics is in marked contrast to its negligi- 
ble effect on the static structure of the melt. The local structure shows 
only a small dependence on the stiffness, and is independent of chain 
length. Calculated correlation functions related to double-quantum 
NMR experiments are in qualitative agreement with experiments on 
entangled melts. A plateau is observed in the dependence of segment 
reorientation on the mean-squared displacement of the corresponding 
chain segments. This plateau confirms, on one hand, the existence of 
reptation dynamics. On the other hand, it shows how the reptation 
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picture has to be adapted if, instead of fully flexible chains, semirigid 
chains are considered. 



1 Introduction 

Modern double-quantum nuclear magnetic resonance (NMR) experi- 
ments aim to understand the microscopic dynamics of polymer seg- 
ments in the melt .111 The polymer dynamics is most often described by 
either the Rouse! or the reptationlla model depending on chain length. 
Both models, however, take the architecture of a specific polymer into 
account only via the so-called Kuhn length. These local features are 
summarized into = j^t, the ratio between the mean squared end- 
to-end distance (R 2 ) and the number of monomers N (multiplied with 
the squared bond length 1%). Statics and dynamics are thus renormal- 
ized onto a Gaussian bead-spring chain of larger but fewer beads of 

(R 2 ) 

the size of the Kuhn length 1% = -jjf-- Experimentally, however, 
there is a qualitative difference between fully flexible polymers like 
poly-(dimethyl-siloxane) (PDMS)El and moderately stiff systems like 
poly- (butadiene) (PB)Q when it comes to the reorientation of chain 
segments. These differences remain after both polymers have been 
appropriately renormalized onto Gaussian chains. The PB system 
shows unusual dynamics which has been taken as an indication of 
local order,H whereas the PDMS data is quite successfully described 
by the standard reptation picture]! From the shape of reorientation 
auto-correlation functions a relatively high degree of residual struc- 
tural order in the presence of entanglements has been deduced for 
PB. The experimental data, however, cannot be interpreted without 
model assumptions. Here, simulations which are validated against 
experimental raw data may contribute to a better understanding.13 

Molecular dynamics simulations are widely anplied to study the 
namics of simplified polymer models in solution^ and the meltH 
To date, only simple models are capable of investigating dynamics of 
long entangled polymer chains in the melt to a satisfactory degree of 
accuracy. Models that allow for atomistic details are limited to shorter 
times and to fewer and shorter chains. 

Recently, we showed that the static local mutual orientation of 
neighboring chains depends on chain stiffness in the amorphous melt, 
wh ereas the chain length has no influence on local static proper- 
tiesOEU In the present contribution, we now extend our investiga- 
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tions to the dynamics of entangled and unentangled melts of poly- 
mer chains with local stiffness. In the following section, the polymer 
model is shortly recapitulated and details of the simulated systems 
are described. In the main part (Section ||), reorientation correlation 
functions are analyzed and compared to theoretical considerations and 
experiments. Section |] relates static chain packing and dynamic chain 
reorientation observables. 



2 Model and Computational Details 

We performed polymer melt simulations using Brownian dynamics 
(BD) of a widely used and well characterized generic polymer modeS§0 
with added local stiffness .Ell All monomers interact via a truncated and 
shifted, therefore purely repulsive, Lennard-Jones potential (Weeks- 
Chandler-Andersen, WCA, potentialEa) 



VwcA(r) = 4e 



<7\ 12 /cr\ 6 1 
r) \r) 4 



, r < V2a . (1) 



Neighbors on the chain are connected by a finitely extensible non- 
linear elastic (FENE) potential which is used for computational effi- 
ciency 

Vfene{t) = -—In^l-^J, ( 2 ) 

r <R = 1.5<7 , a = 30e . 

This yields, together with the WCA-potential, an anharmonic spring. 
Most of our systems have an additional three-body potential to stiffen 
the chain locally 

Vbend = b ( 1 - r% ~ 1,% ri ' l+1 ) , Tij := fj - fi . (3) 



For details of the implementation and parallelization, see ref. 23. 
Throughout this work, reduced units are used with mass m, bead 
diameter cr, and the strength of the WCA potential e set to unity. 
The time unit is t* = &\f^- Temperature is measured in units of e 
by setting Boltzmann's constant ks = 1- The average bond length 
is lb = 0.97o" so that the beads overlap only slightly. 

Systems of up to 1000 chains of length ./V monomers ranging from 
5 to 1000 were simulated at melt density (p* = 0.85<t -3 ) and tem- 
perature T* = 1. The Brownian dynamics algorithm is mainly used 
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to maintain this temperature^ Additionally, it has been shown to be 
very efficient for our modelJiHI The monomer friction coefficient of 
0.5 inverse time units used here means that only processes on a time 
scale well above one of our time units have a meaningful dynamics. 
The persistence lengths l p defined via the decay of the bond direction 
correlation function along the chain backbones varied from one bond 
length up to 5 bond lengths. 

(£(0)n(Z)) = e-'/S (4) 

The distance I is measured along the contour of the chain. The per- 
sistence length is related to the Kuhn length (Section [l]) in the worm- 
like chain model by l p = ^Ir^ In the simulations, l p is controlled 
by choosing appropriate values of b in Eq. ^. Hence, throughout this 
article we refer to systems of different stiffness with their l p rather 
than the corresponding b. The most flexible system (L ~ Z&) has no 
intrinsic stiffness, only the excluded-volume interaction leads to its 
persistence length being non-zero. The persistence length increases 
weakly with chain length due to end effects. Ell The values indicated 
here are the limits for long chains. For N > 10 the persistence length 
is much shorter than the contour length. An overall Gaussian be- 
havior is therefore expected and confirmed by the characteristic ratio 

limTv^oo Re "1f end 6.0 Note that the model would show a nematic 
liquid crystalline phase if the bending stiffness was increased far above 
l p = 5£— I 

All chains of length up to 75 and all flexible chains (l p = 1) relaxed 
fully, as evidenced by the decay of all the Rouse modes. To cut down 
necessary equilibration times, the longer chains (or chains of greater 
stiffness) were initialized as non-reversal random walks whose local 
structure was estimated from simulations of shorter chains: In the 
setup configuration a monomer i and its second neighbor i + 2 are 
not allowed to approach closer than a certain distance. This setup 
procedure reduces theequilibration time substantially while producing 
useful configurations .113 

The end-to-end distance and the gyration radius changed then only 
very slightly in the initial stage of the simulation. Their equilibrium 



values as a function of stiffness were already presented in reference 21 
together with other static observables like structure functions, pair dis- 
tribution functions and local chain orientation correlation functions. 
Some of this data is included in section 0. It is not yet possible to 
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simulate the "slowest" systems (e.g. l p = 5, N = 1000) until the 
final regime of free diffusion is reached. Still, we trust that even this 
system is sufficiently equilibrated for the purpose under study. 

Table g gives an overview of the simulated systems with their re- 
spective Rouse times and simulated times. We are aware that for 
non-flexible polymers the Rouse modes are no longer the true eigen- 
modes (see section ||). Still, the Rouse time is useful as an estimate 
of the relaxation time. For chains of equal length, the Rouse time tr 
would increase linearly with chain extension, i.e. l p , if the friction due 
to neighboring chains were constant. @ However, the relaxation times 
increase even stronger (Table |] and Figure ||). Similarly, for the non- 
flexible chains (l p > 1), the increase of tr with TV 2 as expected from 
the Rouse model is no longer observed. The slowdown is stronger, 
indicating an earlier onset of entanglement influence with increasing 
persistence length (Table [l]). 

The very short chains {N < 15) allow an estimate of the diffusion 
coefficient of chains of a given stiffness, i.e. their mobility in the ab- 
sence of entanglements. The diffusion coefficient does not decrease lin- 
early with chain stiffness as would be expected from the Rouse model; 
for reptation this decrease is quadratic. As entanglement length we 
take the chaiaJength for which the crossover from linear to quadratic 
takes place.c3c3 



3 Chain Reorientation 

3.1 Reorientation Correlation Function 

The main purpose of this work is to investigate the reorientation dy- 
namics of local chain segments in dense melts. This was studied by 
means of the auto-correlation function of the second Legendre poly- 
nomial P2 of chain tangent vectors 

C reor {t) = [3(n(t)n(0)) 2 - l] \ = (p 2 [n(t)u(0)]) (5) 

As chain tangent vectors we take (normalized) vectors connecting 
neighboring beads 

u= (6) 

\n+i - n\ 
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Table 1: Simulated systems: l p : persistence length in monomers; N: chain 
length in monomers; N c : number of chains, t S i m : Simulation times; tr. Rouse 
times (obtained by exponential fit of the lowest 3 to 5 Rouse modes. * partly 
estimated), numbers in brackets for systems where chains are too short; D: 
center of mass diffusion constant (for the systems where free diffusion could 
be reached). Some data of the system with 120 chains of length 350 are also 



presented in ref. |16]. The errors in tr are at most 5% and the errors in QDN 
are about 0.01. 
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Figure 1: Increase of Rouse times with persistence length at N 
dashed line corresponds to an algebraic increase with Z 1 / 54 . 



50. The 



7 



unless noted otherwise. The second, rather than the first, Legendre 
polynomial is taken because its Fourier transform relates directly to 
NMR measurements (Ti experiments) Also the double-quantum ex- 
periments aimed at the chain dynamics are related to this function 
(below). 

3.1.1 Short-time behavior 

For all systems investigated, we have found that the reorientation cor- 
relation function (Eq. ||) consists of two qualitatively different parts. 
For short times, its decay follows a power law (algebraic). At long 
times, the decay is exponential. All characteristics of the reorienta- 
tion correlation function are influenced by the chain length N and 
the chain stiffness l p : the short-time part, the long-time part, the 
time at which the cross-over from algebraic to exponential behavior 
occurs and the function value at this point. We will see that the short- 
time regime is influenced more by stiffness, while the long-time regime 
shows a stronger dependence on chain length. 

The influence of the chain architecture on the short-time behavior 
is illustrated in Fig. [| Reorientation is slowed by increasing the stiff- 
ness at constant length (Fig. ||a) as well as by increasing the length 
at constant stiffness (Figs. ||b and §c for l p = 5 and l p = 1.4, respec- 
tively). The more flexible chains and the shorter chains reach expo- 
nential behavior earlier. At the same time, the short-time algebraic 
process is more "efficient" for the shorter and the more flexible chains, 
i.e. the reorientation correlation function has decreased to a smaller 
value when long-time exponential behavior sets in. These findings can 
be understood if one attributes the short-time behavior to local dy- 
namics and the long-time behavior to the relaxation of larger chain 
segments or, ultimately, to the rotational diffusion of entire chains. 
If the chains are flexible much of the reorientation can be achieved 
by local rearrangement, i.e. without the local tangent vector feeling 
that it is part of a long entangled chain. As stiffness increases, it hin- 
ders the reorientation on local scales, so the reorientation of the local 
tangent vector has to wait for some larger-scale reorientation which 
is exponential. The early crossover to exponential behavior and the 
apparently small efficiency of the short-time process for short chains, 
on the other hand, is due to a faster rotational diffusion of the entire 
chains, as they become shorter. The rotational diffusion begins to 
contribute substantially to the reorientation, before the local process 
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Figure 2: Short-time behavior of time dependent second Legendre polyno- 
mial of next neighbor vectors for different systems, a) N = 200 different 
persistence lengths, b) l p — 5 different chain lengths, c) l p = 1.4 different 
chain lengths, d) System l p — 1.4 for very short times rescaled by jj. 
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can complete. Its time dependence is exponential according to e.g. 
the Debye model. E3 

If the Rouse model were strictly applicable the relaxation time 
of entire chains in the melt (Rouse time tr) should scale with iV 2 . 
Hence, for constant stiffness the reorientation correlation functions 
belonging to different N should coincide if the time axis is transformed 
t — > t/N 2 . Instead, we find empirically that coincidence is achieved 
for t — ► t/N (Fig. ||d). This indicates that even for a small deviation 
from full flexibility (L=1.4) the Rouse model is not appropriate for 
the short-time relaxation, the local process being dominated by effects 
other than connectivity, e.g. stiffness. Additionally the "bead friction" 
enters which also increases sublinearly with stiffness. §1 

In contrast to the local dynamics._there is no chain length influence 
at all on local structural properties .Ej 

As a general result, we note that, while the short-time process is 
influenced by both stiffness and chain length, the influence of stiff- 
ness is much stronger. (This can, for example, be seen by comparing 
Figures ||b and |2|c.) It is, therefore, to be expected that also in ex- 
periments on real polymers the short-time regime will experience the 
influence from the chemical architecture of the polymer, whereas the 
chain length should be secondary. 



3.1.2 Long-time behavior 

The long-time tail of the reorientation correlation function is expo- 
nential to a reasonable approximation (Figure ^). The rate of decay 
depends - as the short-time behavior - on both chain stiffness (Fig. ||a) 
and chain length (Figs. and c). The reorientation time of the ex- 
ponential part with both chain length N and stiffness 
l p (Fig. ||). The dependence on ./V is stronger than in the short-time 
regime, where an approximately linear increase with chain length was 
found (Figure ^d). This is expected because for larger scale processes 
the Rouse model should be appropriate and the overall chain relax- 
ation characterized by the decay of all Rouse modes should scale with 
N 2 . Therefore, the time axis is rescaled in Figures ||b and ||c to high- 
light the remaining deviation from Rouse behavior. As one approaches 
the entanglement regime (the entanglement lengths of the systems are 
compiled in Table [|), an even stronger increase with chain length is 
expected. 

The reorientation time r reor (Figure f|) is much shorter than the 
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Figure 3: Late stage exponential decay of C reor (t) for different systems. Time 
is rescaled by to show differences to Rouse behavior. a)-c) as in Figure § 
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Table 2: Approximate entanglement lengths N e (in number of monomers) 
and respective times r e = Tn(N e ) for the different persistence lengths. These 
were determined from the mean-squared displacements of chains. EE For 
l p = 5 the entanglement time is not defined. 



time for the reorientation of entire long chains. Thus, one may suspect 
that in the case of long entangled chains, it is no longer important 
for the reorientation of a tangent vector, whether the entire chain 
reorients completely. Instead, the relaxation of a shorter part of the 
chain appears to be sufficient. 

This is in line with theoretical expectations: In the limit of in- 
finitely long chains, the local segmental dynamics has to become in- 
dependent of the chain length. The relaxation of a large but finite 
segment of the chain (a few entanglement lengths long) has to give 
enough freedom for the local reorientation to take place. This would 
only be different if the ends of the relevant segment were constrained 
for all times. We were able to corroborate these considerations by 
simulating an entangled melt of fully flexible chains (N = 350, L = 1) 
where the initial algebraic decay of C reor is very fast (» 1800£*)l!iland 
effective in the sense that it reduces C reor (r e ) to less than 0.01 before 
the long-time reorientation sets in. The decay time of the long-time 
process is about 5000t* which is the relaxation time of chain-segments 
of the length of about 60 monomers. This corresponds to about two 
entanglement lengths of the system,0 so we can deduce that, after 
completion of short-time relaxation, only relaxations on a length scale 
of up to the order of the entanglement length are relevant. 

3.1.3 Reorientation of medium-size chain segments 

It is also of interest to analyze the reorientation of longer chain seg- 
ments. As the orientation vector of a segment of length d we take the 
unit vector between two monomers whose indices on the chain differ 



12 



10 c 



icr 



10 



o 

CD 



10' 



10' 




10' 



10 



100 

N 



1000 



Figure 4: Reorientation times r reor of local segments as a function of chain 
length and stiffness. These are calculated from exponential fits to the long- 
time part of the reorientation correlation function Eq. [|. The dashed line 
indicates an increase with iV 2 ' 3 for l p — 5 whereas the solid line indicates 
N 2 - 9 for L = 1.4. 
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Figure 5: Reorientation of chain segments of different length d for the system 
with length N = 200 and persistence length l p = 5: d—1, 3, 5, 9, 11, 13, 19, 
39, 79, 119 and 199 (end-end vector) from bottom to top. 



by d 



-, r i+d r i /_-, 

u d = -p =rr . (7) 

\ r i+d ~ r i\ 



Thus d = 1 represent the vectors connecting neighbors discussed up 
to now. Using simple topological arguments, the dependence of the 
exponential part of the reorientation correlation function (Eq. ||) of u d 
on d has been predicted.^ The value of the reorientation correlation 
function at any given time is proportional to (3d/^)/(5A r e /fe), where 
iV is the chain length, N e the entanglement monomer number, and 
Ik the Kuhn length. This relation holds for Ik < Nit, < N e lb. As a 
measure of the amplitude of the exponential part of the reorientation 
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Table 3: Fitted exponential decays to the curves in Figure [5| (N — 200, l p = 
5) in the time domain between t = 100000 and 300000. The amplitude (3 is 
defined by the following relation C reor (t) — > (3e~ t / Treor , t — > oo, i.e. it is the 
fictitious intersection of the exponential long-time curve with the y-axis. 



correlation function we define the ordinate intercept (3 obtained by 
fitting an exponential to the long-time tail of the reorientation corre- 
lation function and extrapolating back to t = 0. If the assumptions of 
the theory were true, (3 would be proportional to d. In Table ||], how- 
ever, (3 is seen to increase monotonically but sublinearly with d. The 
data of Table || correspond to a non- flexible chain (l p = 5). Hence, 
it is obvious that the arguments of topological entanglement are not 
sufficient to explain the behavior of stiff chains. The reason for this 
is most likely that for the l p = 5 system Ik and N e lf, are of the same 
order of magnitude, as shown elsewhere,@ so that the above condition 
is not fulfilled. 

The long-time exponential parts of the reorientation correlation 
functions belonging to different d are almost parallel (Fig. |5|). They 
have relaxation times T reor of similar magnitude, although the T reor 
seem to increase slowly but monotonically with d (Table |3|). 

3.2 Comparison to Double- Quantum NMR Ex- 
periments 

The direct experimental observable in double-quantum (DQ) experi- 
ments like the ones performed by Graf et alB is 

C DQ (t) = (p 2 [Bm]P2[Bu(t)]) ■ (8) 
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The vector u is a unit vector along an atom-atom distance vector 
which is usually not parallel to the chain tangent vector. If enough 
of these are available, then the Cdq of the backbone can be recalcu- 
lated .§1 B is a unit vector parallel to the external magnetic field in 
the NMR experiment. We can choose B = e z for convenience because 
amorphous melts are rotationally invariant. Cdq is proportional to 
C reor if the vectors u are isotropically distributed with respect to the 
field 

C DQ (t) = \c reor {t) . (9) 
o 

can be used for the comparison to experiments. How- 
ever, absolute values cannot be compared, since the experimentally 
detectable C reor (0) is reduced from 1 to a value S by very fast mo- 
tions of internal degrees of freedom not present in our bead-spring 
model. For the C=C double bond in polybutadiene, S is found to 
be 0.24, for exampleB For comparison, we therefore normalize both 
curves to C reor (0) = 1. 

Ball et al. have derived expressions for the DQ correlation func- 
tions assuming the reptation model and infinitely long chains.il They 
predict an algebraic decay with different exponents in the different 
dynamic regimes of the standard tube model. In the time interval 
between the entanglement time r e and the Rouse time tr, for which 
the inner degrees of freedom of the chain are relaxed, a i -1//4 regime 
of C reor is expected. Later, in the regime where the chain as a whole 
reptates in its tube, a t _1//2 behavior should be found. The exponent 
of C reor should be the negative of that of the monomer mean-square- 
displacement in the same dynamic regime. Algebraic fits of the short- 
time part of the decay curves (linear region of the double logarithmic 
plot cf. Fig. §a-c) yield exponents k shown in Table ||. The systems 
under study are not very long compared to the infinite chain limit as 
they are at most about 30 times the entanglement length (except for 
N = 1000, l p = 5). Our exponents k come therefore closer to the t^ 1 / 2 
dependence. The system with persistence length l p = 5 is the most 
strongly entangled.cj It is found to reorient slowest with an exponent 
between 0.25 < k < 0.5. The system with persistence length l p = 1.4 
shows an algebraic decay faster than t~ 1 / 2 . This is probably because 
it is so weakly entangled, that the effects of entanglements are just 
starting to play a role. The exponents decrease systematically with 
persistence length, which, as discussed earlier, indicates an increasing 
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Figure 6: Comparison of C reor from our simulations (l p = 5) with the exper- 
iments on polybutadiene of ref. Both are scaled to C(0) = 1. The time 
axis for the simulation data is scaled empirically as an independent mapping 
to experimental times is not possible. The experimental times are measured 
in entanglement times which are derived from viscosity measurements. 



17 



lp 


K 


1.4 


0.68 


3.0 


0.50 


5.0 


0.40 


5.0* 


0.29 



Table 4: Algebraic fits (t~ K ) of the decay of double-quantum 

correlation functions Cdq for iV = 200 (see text and Figure H). 
* The bottom line for l p — 5 has chain length TV = 1000 but the system 
is not equilibrated. 



degree of entanglement. The dependence on the degree of entangle- 
ment is supported by the very low exponent (re = 0.29) for a system 
with l p = 5 and chain length N = 1000. 

A i~ 1//4 dependence of the correlation was found experimentally 
by Graf et al.u for a system with very long and therefore highly en- 
tangled PB chains (76 entanglement molecular weight M e ). They 
also observed a power-law with t -1 / 2 for a system with a molecular 
weight of HM e . Figure || compares directly simulation (at l p = 5) 
and experiment; the ratio between the lengths in simulation and ex- 
periment is not exactly the same, it is, however, only important to be 
slightly or far above M e . The agreement shows that the simulations 
reproduce well the exponents found in experiments. To achieve this 
agreement, the time axis of the simulated correlation functions has 
been rescaled empirically by 0.153 and 0.5 for TV = 50 ~ 8N e and 
N = 1000 ~ 160iV e respectively. This scaling may be used to infer a 
mapping to experimental times. 

3.3 Interdependence of reorientation and trans- 
lation of segments 

If the reptation model holds the reorientation process is coupled to the 
translation of the polymer in its tube. A useful relation to monitor is, 
therefore, the reorientation correlation function of the chain tangent 
vector versus the mean-squared displacement of the monomers defin- 
ing it, irrespective of the time. This relation has to be averaged over 
a finite time window 2t av , which is centered at some time t m and does 
not necessarily start at t = 0. 
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C reor {Ar 2 ) = (P 2 [u(0)u(t)]) , (10) 

Ar 2 = /[r(t)-r(0)} 2 \ (11) 

\ / t m t av ■ l t rn -\-t av 

Both C reor and Ar 2 depend parametrically on t. If, during t, the 
tube relaxes (reorients) then C reor is zero. Any deviation from zero 
indicates that the orientation is correlated over the time interval cor- 
responding to the displacement. In our analysis we have ruled out 
possible artifacts from the translation of the system as a whole which 
could be present in Brownian dynamics. 

In Figure |7|a, it is seen that this function at short t (0 < t < 30000) 
does not decay to zero, but shows a plateau, whose value depends on 
the stiffness: stiffer chains have a higher residual correlation. The 
presence of a plateau is a consequence of the finite length of the chain: 
Every finite polymer has a trivial residual static orientation correla- 
tion between distant chain segments in the direction of the end-to-end 
vector. If the motion of the chain is predominantly along the fixed 
tube, this residual correlation is also visible in the dynamic C reor (Ar 2 ) 
shown here, since one chain segment samples the very same fixed tube 
at different times. We have seen in the preceding sections that stiffer 
chains have a higher reptation component in their dynamics. Hence, 
it is no surprise that they exhibit a larger residual correlation. 

Figure |^b shows, for the most interesting case l p = 5, how the 
C reor (Ar 2 ) depends on the position of the time window t m . With 
increasing t m , the reorientation correlation at Ar = goes to zero. 
At the same time, intensity moves to larger Ar 2 , so that eventually 
a maximum at Ar 2 > develops. This observation is explained by a 
scenario that includes, in addition to reptation, a diffusive or rather 
subdiffusive translation of entire chain segments through space, with- 
out significant reorientation of these segments: When a chain segment 
has reptated along the tube and comes back to its former part of the 
tube and hence its former orientation, it finds that this part of the 
tube itself has translated in the meantime. On the other hand, when 
it returns to its former absolute position (Ar 2 = 0) it is now in a com- 
pletely different part of the tube and has a different orientation. This 
picture of short-scale transverse translation of stiff tube s e g ments has 
been borne out by visualizations of individual chains£3E3 The time 
dependence of the plateau value contains information about the sta- 
bility of the initial neighborhood. It measures the "similarity" i.e. 
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Figure 7: Reorientation correlation function depending on the mean-square 
displacement of monomers (N = 200) a) Different systems at < t < 30000. 
A running average was applied after the initial decay in order to show the 
possible plateau more clearly on all the curves in Figure a) and for the curve 
of t m = 75000 in Figure b) So the apparent final increase results from low 
statistics and is not meaningful, b) l p — 5 three different time intervals 
t m and comparison to static correlation. For reasons of statistics we always 
average over 2t av = 30000, i.e. t m = 15000 means < t m < 30000. The 
static case is measured against topological distance, because distance in real 
space is not appropriate as the chain can fold back which leads for packing 
reasons (Section |j) to perpendicular alignments. 
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correlation of the neighborhoods at different points in time as the 
chain returns. The neighborhood changes with time on a much longer 
time scale. 

These results support the presence of reptation in our systems, as 
the chains come back to their former surrounding which has undergone 
only small changes in the meantime. As this memory effect preserves 
information about orientations, a tube picture is a suitable concept. 
However, the chains do not behave simply as the standard reptation 
picture would suggest. The reptation is considerably modified by their 
stiffness. Stiffer chains reptate in a more pronounced way, i.e. they 
follow the primitive path of the tube more closely as the stiffness 
suppresses the transversal motions efficiently This leads to a higher 
degree of orientation memory for chains of the same length (Figure 0a) . 

4 Connection to Structure 

The preceding section has shown that both the stiffness l p and the 
length TV have an effect on the dynamics of polymers alreadv on the 
local level. In this section, we briefly review earlier resultstH about 
the local structure in polymer melts and how it is influenced by the 
chain architecture. This is done not only for comparisons within the 
model system. NMR experiments on melts have so far only been 
able to study the local dynamics. Any information on the structure 
had to be deduced from the dynamics using models. In contrast, 
the simulations of this work can be analyzed independently for both 
structure and dynamics. We, therefore, have an example case for 
which the assumptions can be checked that are used to analyze NMR 
experiments. 

Of particular interest has been the question whether or not neigh- 
boring polymer chains are in any way aligned.E3 We therefore concen- 
trate here on the static orientation correlation function OCF 

OCF(r) = (p 2 [u l u 2 )(r)Y (12) 

This function measures the mutual orientation of tangent vectors (de- 
fined in Eq. |6|) of segments belonging to two different chains 1 and 2 
as a function of their distance r. The second Legendre polynomial is 
used again, this time because our polymer chains have no direction, i.e. 
head and tail are equivalent. The OCF is 1 for parallel orientation, 
— 1/2 for perpendicular orientation, and for random orientation. 
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Figure 8: Comparison of the static orientation correlation of a) chains of 
different lengths (fully flexible system, l p — 1) and b) chains of different 
stiffness (N = 50) 

The detailed discussion of the various OCFs is given elsewhere. EE 
Here we only note that the OCF is a strictly local property. The chain 
length TV has no influence whatsoever on the short-range mutual ori- 
entation of two chains, even if one N is below the entanglement length 
iV e ~ 32 and the other above (Figure ||a). The influence of the stiff- 
ness l p on the structure is clearly visible (Figure ||b), but small. We 
may conclude that the two parameters ./V and l p , which both influ- 
ence significantly the local reorientation dynamics, have little (l p ) or 
no (N) effect on the local packing of chains. The fact that the chain 
length is not important for the local structure means that entangle- 
ments cannot be important either. This is yet another manifestation 
of the entanglement length N e being a purely dynamical quantity. 

Local packing is, however, strongly influenced by another local 
quantity (which, in contrast, contributes little to the dynamics^) , 
namely the excluded volume of the monomers. An example of this is 
shown in Figure |. Here, the OCF for a flexible chain (l p = 1, N = 50) 
is overlaid by the OCF for dimers on a randomly perturbed lattice. 
The monomers occupy fee lattice sites whose positions were randomly 
displaced by small amounts to emulate finite temperature. The OCF 
was then evaluated between all possible pairs of dimers which do not 
share a common atom. Although the OCF for the dimers is much 
more accentuated than for the amorphous melt one can clearly see 
the short-range orientation of dimers shining through in the OCF of 
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Figure 9: Comparison of the static orientation correlation of simple bead- 
spring chains against orientation correlation function of dimers on a per- 
turbed fee lattice, the dotted line is a spline curve as a guide to the eye. 
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the melt. 

One may, therefore, conclude that local structure and local dy- 
namics in polymer melts are dominated by different properties of the 
polymer. For this reason, it may be difficult to infer one from experi- 
mental results on the other .0 



5 Conclusions 

The reorientation of short segments in polymer chains in dense melts 
is governed by two subsequent processes. The fast one leads to an 
algebraic decay of the reorientation correlation function and the slow 
one to an exponential decay. The correlation times of both depend 
on the chain length as well as on the chain stiffness. Increasing chain 
stiffness leads to a strong slowing of the reorientation on both time 
scales. 

A qualitative comparison of our reorientation correlation function 
showed that the power laws of the reorientation correlation function 
measured in double-quantum NMR experiments of systems not too 
far above the entanglement molecular weight could be reproduced. 
Therefore, our simplistic model, which is probably the simplest possi- 
ble to incorporate stiffness and excluded volume, successfully describes 
the qualitative features of the dynamics. With our results, thus val- 
idated against experimental data, the reorientation correlation func- 
tions C reor {t) can be regarded as meaningful. In contrast, our simple 
model cannot explain two other experimental observations. The initial 
plateau of C reor (t), the so-called dynamical order parameter,! of the 
experiment is not found, and in experiment the difference between per- 
sistence length and entanglement length is larger. Detailed atomistic 
models are probably necessary, in order to capture these features 

One main result of this investigation is that, for the local reorien- 
tation to relax slowly, the chains have to be both stiff and long. Local 
stiffness ensures that the memory of orientation is not already lost 
during the fast algebraic process. On the other hand, the chains have 
to be above the entanglement length for the slow exponential process 
to extend into the experimentally observable regime (milliseconds in 
NMR). Although in principle both processes are always present, they 
have to span a big enough range in intensity and time, in order to be 
detectable in experiments. In agreement with our results, it has been 
found previously that the chain length has only little influence on the 
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local reorientation, as long as it is below the entanglement length. In 
united atom polyethylene below the entanglement length local orien- 
tation relaxes like a stretched exponential, whereas the relaxation of 
the end-to-end vector is well described by the Rouse model.Ej 

Although the dynamic reorientation correlation functions are found 
to be in qualitative agreement with experiments and theoretical pre- 
dictions, our work shows that this does not necessarily imply an in- 
crease of static order due to topological entanglements of the poly- 
mer chains. The static local order increases with chain stiffness but 
does not at all depend on chain length. The entanglement length 
has emerged as a central quantity to explain the dynamics of flexi- 
ble and especially stiff chains. It is not easily determined uniquely as 
there are several different definitions which lead to different values. 
This difficulty becomes even worse for chains with intrinsic stiffness. 
Nonetheless, we can say that the entanglement length in any definition 
decreases dramatically as the persistence length is increased. Evidence 
for this is found in the chain length dependence of the center of mass 
diffusion coefficient and in the segment size governing the long time 
relaxation of local order. In the system with l p = 5, it is question- 
able whether any length scale can be described by Rouse dynamics. 
The very local scales up to the persistence are dominated by bending 
modes and the very big length scales are dominated by entanglements. 
As NJb approaches l p , the Rouse regime disappears between these two 
extremes. 

A renormalization of the local scale properties onto an effective 
monomer or Kuhn segment is possible for static structural aspects as 
there is only one relevant length scale, namely the persistence length. 
However, this renormalization fails for dynamical properties of stiff 
polymers because two length scales (and the associated time scales) 
interdepend.^! In dynamics one encounters, in addition to the persis- 
tence length, the entanglement length which describes the topological 
constraints imposed by non-crossability of the chains. When the two 
scales come into the same order of magnitude, new behaviors emerge 
which can not be deduced by renormalizing to Rouse dynamics or 
other simple models. The fact that the dynamics is not described ap- 
propriately by analytical theories implies also that the connection of 
translational and rotational dynamics is not a priori known. This has 
to be kept in mind by interpreting NMR experiments which predom- 
inantly measure reorientation. 

The concept of reptation is supported by the existence of a time- 
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dependent plateau value of reorientation as a function of the length of 
the diffused path in our simulations. Reptation is more pronounced 
if the chains are stiffer because of both the intrinsic stiffness and the 
stronger entanglement. They lead to the chain being more closely 
confined to the primitive path of the tube. 

An analysis of the dependence of N e on l p and further consequences 
for chain and monomer motions will be discussed in more detail else- 
where. H 
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A List of Symbols 



a 


strength of FENE potential 


B 


magnetic field axis 


b 


strength of bending potential 


P 


amplitude of reorientation correlation function 


Cdq 


double quantum correlation function 


C re or 


reorientation correlation function 


D 


center-of-mass diffusion constant 


d 


chain segment length 


e 


nonbonded interaction strength 


k B 


Boltzmann's constant 


K 


decay exponent for correlation function 


I 


distance along chain contour 


h 


bond length 




persistence length 


h< 


Kuhn segment length 


m 


monomer mass 


M e 


entanglement molecular weight 


N 


chain length in monomers 


N c 


number of chains 
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entanglement monomer number 


OCF 


orientation correlation function 


P2 


second Legendre polynomial 


R 


characteristic length of FENE potential 


-^ctid, — end 


end-to-end distance 


It/10 i>y* 


gyration radius 


r 


distance 


r 


position vector 




r 


density 


s 


dynamical order parameter 





monomer diameter 


T 


temperature 


1 


spin-lattice relaxation time 


t 


time 


t* 


time unit 


t a V 


averaging time 


tm 


point of time for correlation of orientation and translation 


tsim 


simulation time 


T e 


entanglement time 


TR 


Rouse time 


Treor 


reorientation correlation time 


U 


tangent unit vector 


u d 


tangent unit vector for segment length d 


V 


interaction potential 
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